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Matrix Multiplication 

Complete Program 

198 - 159 + (N12) 

Up to 12 places . This will depend on the number of places 

to be punched. 

Let A he an (i x j) matrix and B he an (j xk) matrix. The 

number of elements of these matrices must be such that 

j k< 10,240 and j < hl2 . i is arbitrary. 

Let JZ , m, and n be the number of decimal place in the A, B, 

and product matrices respectively (k jkm + 18 jk -f 80 k + 500) 

millisec. to input B, and (3.3 ijk + b.lj£ + 6.6 ikn + 200 ik) 

millisec. to input A and perform calculation with output. 

0, 1, 2, 3, h, 5. 

This program will take an (i x j ) matrix A and multiply it 

by a (j x k) matrix B yielding an (i x k) matrix C. No 

scaling is necessary other than that each element of the 

matrices be less than one. If overflow occurs, the program 

records it and punches out the proper element of matrix C 

accordingly. Each element of C which is larger than one is 

punched as a signed integral part, a decimal point, and an n 

place fractional part. Otherwise the element of C is 

punched as an n place signed fraction. The elements of the 

matrix C are preceded by three constants i, k, and n, where 

i indicates the number of rows of C5 k, the number of 

columns, and n, the number of places printed. Each row 

is terminated by the letter N. Thus if all the elements of 

C are less than one, the matrix C can be read back into the 

computer again. 
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METHOD OF OPERATION 



PREPARATION OF MATRIX TAPE 

Prepare the matrix A tape by first punching the constants 
i, j, and n (see note l) where i denotes the number of 
rows of A| Oi "the number of columns of A; and n, the number 
of decimal digits of c., punched. These constants are to 
be separated from each other by a fifth hole character 
(e.g. space). After these constants, punch the elements 
of A in a rowwise fashion. Each row of A is terminated 
with the symbol N (see note 2)* Punch the matrix B in a 
similiar fashion (see note 3}» 

After the program has been read, a 2k stop occurs. Place 
the B matrix on the reader. Upon bypassing this step the 
matrix B is read and placed on the drum (see note k)* A 
30 stop follows . Now place matrix A in the reader and 
bypass the stop. The matrix A is read a row at a time and 
k elements of matrix C are punched for each row of A (see 
note 5)» When the matrix multiplication has been completed, 
the program will stop on a 30 order. If one bypasses this 
stop with the black switch, a new matrix A will be read and 
multiplied with the old matrix B. If one bypasses this 
stop by means of the white switch, both a new matrix B and 
a new matrix A will be read and multiplied together (see 
note 6 ) . 

1. The constants i and j need not coincide with the actual 
number of rows and columns of the matrix but may be less 
If this is the case, only the upper left hand (ixj) 
submatrix is used. 

2. If instead of using an N as a row terminating symbol, 
one used an J, then the program coming upon it will 
stop. When this stop is bypassed by the black switch, 
the program will resume input of the matrix. This 
enables one to have the matrix on several tapes. 

3* The last printing constant n to be read governs the 
number of places printed. The previous printing 
constants have no effect on the printing of the current 
matrix . 
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k. If the FF 03N stop occurs, this signifies that 
SADOI has played back the program improperly. 

5. If the FF 03J stop occurs, this signifies that 

the matrix B has been read off the drum incorrectly. 
This stop will occur if the constant j preceding 
the matrix A does not agree with the constant j 
preceding the matrix. B. 

6. Automatic operation may be achieved by placing 
the black switch in the ignore position. However, 
certain additional characters must be punched between 
matrices, i.e. after the last N of some matrix, yet 
always prior to the three constants on the following 
matrix. These characters have the following meanings j 
N or J; Causes the program to stop on an OF 

instruction. By working the white switch 
up and down, one returns to the point in 
the program which precedes the reading of 
this character. 
F or Lt Causes the next matrix input to be treated 
as a B matrix, rather than a new A matrix 
as would be the case if no F or L separates 
the two matrices . 

Thus, if one arranges his matrices and inserts the 
character F or L properly one may do a sequence of 
A*B multiplications with just new A matrices or with 
both new A matrices and new B matrices without stopping. 
If now an N or J character follows the F or L character, 
an OF stop will occur. -3/6/61 
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